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We study the dynamics of ordering in ferromagnets via Monte Carlo simulations of the Ising 
model, employing the Glauber spin-flip mechanism, in space dimensions d = 2 and 3. Results for the 
persistence probability and the domain growth are discussed for quenches to various temperatures 
(T/) below the critical one (Tc), from different initial temperatures Ti > T^. In long time limit, for 
Ti > Tc, the persistence probability exhibits power-law decay with exponents 9 ~ 0.22 and ~ 0.18 
in d = 2 and 3, respectively. For finite Ti, the early time behavior is a different power-law whose 
life-time diverges and exponent decreases as Ti ^ T^. The crossover length between the two steps 
diverges as the equilibrium correlation length. Ti = Tc is expected to provide a new universality 
class for which we obtain 6 ~ 0.035 in d = 2 and ~ 0.10 in d = 3. The time dependence of the 
average domain size I, however, is observed to be rather insensitive to the choice of Ti. 


I INTRODUCTION 


When a homogeneous system is quenched below 
the critical point, the system becomes unstable to 
fluctuations and approaches towards the new equilib¬ 
rium via the formation and growth of particle rich 
and particle poor domains [ll-y. In such nonequi¬ 
librium evolutions, over several decades, aspects that 
received significant attention are the domain pattern 
rate of domain growth [p 1^ 15j| , persistence 
Sfllilil, nil and aging [23 - E3| . Avera^size, £, of 
domains typically grows with time (t) as (1] 


r. 


( 1 ) 


The value of the exponent a for nonconserved order- 
parameter dynamics sm, e.g., during ordering in an 
uniaxial ferromagnet, is 1/2, in space dimension d = 2. 
In addition to the interesting structures exhibited by 
the domains of like spins (or atomic magnets) in a 
ferromagnet, the unaffected or persistent spins also 


form beautiful fractal patterns 


. Typically, 


fraction of such spins, henceforth will be referred to 
as the persistent probability, P, decays as 


t 


-e 


( 2 ) 


with [20|, l2l| 0 having a value ~ 0.22 for the Ising 


model in space dimension d = 2 and ~ 0.18 in d = 3. 

The values of the exponents mentioned above are 
understood to be true for the perfectly random initial 
configurations, mimicking the paramagnetic phase at 
temperature T = oo. Another relevant situation is to 
quench a system from finite initial temperature (Ti) 
with a large enough equilibrium correlation length 
However, this problem has received only occasional at¬ 
tention though experimentally very relevant. 

In this context, the behavior of the two-time equal 
point correlation function, relevant in the aging phe¬ 
nomena, was studied 3l|, l32| in d = 2 for Ti = T^, 
the critical temperature. It was pointed out that such 


quenches would form a new universality elass and was 
shown that the decay of the above correlation was sig¬ 
nificantly slower for Ti = T^ than Ti = oo. In view 
of that a slower decay of P is also expected. Apriori, 
however, it is unclear whether there will be quantita¬ 
tive similarity between the degree of changes in the 
two quantities. On the other hand, the behavior of 
P and £ are expected to be disconnected 3^. Nev¬ 
ertheless, the rate of growth of £ may be different for 
Ti = Tc and = oo, at least during the transient pe¬ 
riod. In this paper, we address the Ti dependence for 
persistence and domain growth in a ferromagnet, via 
Monte Carlo (MC) simulations [s^ of nearest neigh¬ 
bor Ising model [3^ 


H = -J^ S^Sf, Si = ±l, J >0, (3) 

<ij> 


in d = 2 and d = 3, on square and simple cubic lat¬ 
tices, respectively. 

Starting from a high value, as Ti approaches T^ 
2.27J/kB in d = 2 or 4.51 J/A:b in d = 3, ks being the 
Boltzmann constant], a two-step decay in P becomes 
prominent, except for Ti = Tc. For the latter, power- 
law behavior with exponents much smaller than the 
ones observed for quenches from Ti = oo lives for¬ 
ever. In addition to identifying these facts, a primary 
objective of the paper is to accurately quantify these 
decays. For the domain growth, on the other hand, 
we do not observe a modification to time dependence 
with the variation of Ti, almost from the very begin¬ 
ning. 

The rest of the paper is organized as follows. In 
section II we briefly describe the methods. Results 
from both the dimensions are presented in section III. 
Section IV concludes the paper with a summary and 
outlook. 
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II METHODS 


The nonconserved order-parameter dynamics in 
the MC simulations have been incorporated via the 
Glauber spin-flip mechanism [^ . In this method, 
a randomly chosen spin is tried for a change in sign 
which is accepted according to the standard Metropo¬ 
lis algorithm [^. We apply periodic boundary con¬ 
ditions in all directions. Time is expressed in units 
of MC steps (MCS), each MCS consisting of L'^ tri¬ 
als, L being the linear dimension of a square box. We 
have computed (. from the domain size distribution, 
Pd{td,t), as [TSI 


^{t)= idPdild,t)d£d, 


( 4 ) 


where id is calculated as the distance between two suc¬ 
cessive interfaces in any direction. All lengths are ex¬ 
pressed in units of the lattice constant a. We present 
the results after averaging over multiple initial config¬ 
urations. This number ranges from 20 (for L = 1024) 
to 200 (for L = 400) in d = 2 and 10 (for L = 256) to 
50 (for L = 64) in d = 3. The initial configurations 
for Ti close to Tc were carefully prepared via very long 
runs. At Tc, depending upon the system size, length 
of such runs varied between 5 x 10® to 10® MCS. 



FIG. 1. Upper panels show snapshots during the evolu¬ 
tion of the Glauber Ising model with Ti — oo, Tf = 0 
and L = 512. The black regions represent domains of up 
spins. The lower panels show the unaffected spins, marked 
in black, corresponding to the evolution snapshots above 
them. These results correspond to d = 2. 


Ill RESULTS 

A. d = 2 

Growth of the domains have been demonstrated in 
the upper frames of Fig. [T]for the system size L = 512 
in d = 2. In this figure we show snapshots from two 
different times during the evolution of the Glauber 
Ising model. In the lower frames of the figure, we show 
pictures marking only the persistent spins. Beautiful 
patterns are visible. These results correspond to a 
quench from = oo to Tf = 0. 

Plots of P, for Ti = oo and few different values of 
Tf are shown vs t in Fig. [5] The data for T/ = 0 
and 0.25rc are consistent with each other and follow 
power-law, the exponent being 0 ~ 0.22. The flat 
behavior at the end is due to the finite-size effects. 
This value of 9 is consistent with the observation by 
Manoj and Ray [2^. However, for hig her values of 
Tf, as also previously observed [l^ ll9|, the decay is 
not of power-law type. This is thought to be due to 
thermal fluctuation. When this fluctuation is taken 
care of, we observe 9 ~ 0.22 for all the values of Tf 
included in Fig. [21 in agreement with Ref. [21]. It 
is thought that persistence and domain growth are 
not strongly connected to each other. Interestingly, 
different behavior in Fig. |2]for T > 0.25Tc is strongly 
reflected in the domain growth also. For brevity we 
do not present it here. 



FIG. 2. Plots of persistence probability, P, vs, time, on a 
log-log scale, for quenches from T = oo, with L = 512, in 
d = 2. Four different values of T/ are included. The solid 
line there has a power-law decay with exponent 0.22. 

In Fig. IHKa) we show P vs t plots, on a log-log 
scale, for quenches to Tf = 0, from a few different 
values of Ti, all for the same system size P = 512. It 
appears that, in the long time limit, for Ti > Tc, the 
decay is power-law, with the same exponent 9 ~ 0.22. 
However, crossover to this functional form gets de¬ 
layed as Ti approaches T^. In fact, in the pre-crossover 
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FIG. 3. Log-log plots of P vs t, for quenches from different 
values of Ti{^ Tc), to Tf — 0. Continuous lines there 
correspond to power-law decays with exponents 0.22 and 
0.04. (b) Instantaneous exponents 6i are plotted vs l/£, 
for the quenches in (a), excluding Tt = Tc case. Here we 
have included only the late time behavior. The continuous 
lines in this figure are linear fits, providing the values of 6 
as the ordinate intercepts. All results are from simulations 
in d = 2. 


regime, another power-law decay, with smaller value 
of 9, becomes prominent with the decrease of Ti. Such 
a slower decay becomes ever-lived for Ti = Tc. 

In Fig. [3jb), we present the instantaneous expo¬ 
nent, 9i, calculated as |l4 Il5| 


9 ,. = - 


dlnP 

dint ’ 


(5) 


vs l/£, with the objective of accurate quantification 
of the second step of the decays for Ti close to, but 
greater than, T^- For the abscissa variable we have 
adopted l/£, instead of 1/t, to visualize the long time 




FIG. 4. (a) Instantaneous exponents 6i are plotted vs l/£ 
for two of the quenches in Fig. El a). Here we have focussed 
on the first step of the decays, exponents for which are 
obtained from the flat regions, marked by the horizontal 
solid lines. These values of 9 are plotted, vs e = Ti — Tc, 
in the inset of the figure. The continuous line there is 
a power-law fit (see text), (b) Plots of P vs t, for two 
different system sizes, with T = Tc and Tf — 0. The 
upper inset shows 9i vs 1/i for three different values of L. 
The dashed line in this inset is a linear extrapolation using 
data in the small £ region. The L-dependent exponents 
obtained from flat regions of the plots (see the horizontal 
solid lines) are plotted vs 1/L in the lower inset. The solid 
line there is a linear fit. All results correspond to d = 2. 


limit better. Data sets for each Ti are fitted to linear 
functions to obtain the exponents. Within statistical 
error, for all the presented temperatures, it appears 
that 9 is consistent with the quench from T^ = oo to 
Tf = 0. Nonzero slopes in these plots are due to the 
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presence of off-sets at the crossover times. 

Next we move to identify the exponent for Ti = 
and Tf = 0. In Fig. [HJa), it appears that the 
Ti = Tc data are reasonably consistent with 6 = 0.04. 
Nevertheless, before the final finite-size effects appear 
(showing flat nature at very late time), there is a faster 
decay, albeit for a brief period. This can well be due 
to the fact that for a finite system, ^ is not infinite 
at T = Tc, effectively implying that the initial con¬ 
figurations are prepared away from Tc- Thus, in this 
problem finite-size effects have two sources. One com¬ 
ing from the finiteness of the equilibrium correlation 
length, other being faced when the nonequilibrium do¬ 
main size is close to the system size. Thus, a quan¬ 
tification of the exponent d, for Ti = Tc, via finite-size 
scaling 37|, becomes a challenging task. However, we 
appropriately take care of the shortcoming below, in 
various different ways which provide results consistent 
with each other. 

In Fig. lUJa) we show the instantaneous exponents 
6i, vs l/£, with the objective of quantifying the first 
step of the decay, for two values of Ti, close enough to 
Tc- As demonstrated, from the flat regions we identify 
the exponents for the early time or the first step of 
the decay for various values of Ti. These numbers 
are plotted in the inset of this figure as a function of 
e = Ti — Tc- The continuous line there is a fit to the 
form 


e{T,)=0{Tc)+Ae^, ( 6 ) 

providing 6{Ti — Tc) = 0.034, A = 0.15 and x = 0.54. 
Recall that this value of 6 is the only decay exponent 
for T, = Tc. 

To verify the above value of 6 further, in Fig. lUJb) 
we present an exercise with different system sizes. In 
the main frame of this figure, we present P vs t data, 
for Ti = Tc, from two different values of L. It is seen 
that with the increase of the system size, there is a 
tendency of the data to settle down to a power-law for 
a longer period of time, following a marginally faster 
decay at early time. In the upper inset of this figure 
we show 6i vs l/£ for three different system sizes with 
Ti = Tc. The early time behavior appears linear, ex¬ 
trapolation of which leads to 0 ~ 0.029. However, 
if the data in the main frame is closely examined, 
as already mentioned above, this part corresponds 
to the preasymptotic behavior, thus, should be dis¬ 
carded from the analysis. Actual asymptotic expo¬ 
nents should be extracted from the flat regions of the 
plots. The numbers obtained from these flat parts, 
as discussed, differs due to the finite-size effects and 
thus, should be extrapolated to L = oo appropriately. 
These values are plotted in the lower inset as a func¬ 
tion of 1/L. A very reasonable linear fit (see the solid 
line) is obtained, providing 0 = 0.035. On the other 
hand, a nonlinear fit appears to be nearly quadratic 


providing 6 = 0.037. From all these exercises we con¬ 
clude that 0(Ti = Tc) = 0.035 ± 0.005. This picture 
remains true for quenches from Tc to nonzero values 
oi Tf, if thermal fluctuation effects are appropriately 
taken care of. On this issue of thermal fluctuation, 
here, as well as for Ti = oo, our studies were restricted 
to only small system sizes. 

The decay of the two-time correlation is also of 
power-law type. For quenches from Ti — Tc, the value 
of the exponent gets reduced by a factor ~ 10, com¬ 
pared to Ti = oo. In the present case the reduction 
factor is ~ 6.3. While there may be connection be¬ 
tween the two phenomena, but a search for matching 
between the two factors may not be justified. This 
is because, only spins with odd number of flips are 
important for the decay of the two-time correlation 
function. 

It is certainly important to ask, if, like the decay 
of the persistence probability and the two-time cor¬ 
relation [^, the growth of the average domain size 
also exhibits initial temperature dependence, at least 
at the transient level. Here we only examine the cases 
Ti = oo and Ti = Tc, for quenches to Tf = 0. In Fig. 
[5] we present the i vs t plots for these two cases using 
a log-log scale. Both the data sets appear to grow 
slower than This can well be due to the presence 
of significantly big initial length Iq, which we examine 
below. While from this figure it is difficult to identify 
any difference in the growth exponent between the two 
cases, there certainly exists difference in the finite-size 
effects, noting that L = 512 in both the cases. 

To learn better about the exponents, in the inset 
we present the instantaneous exponents 14l. Il5| 


d\n£ 
dint ’ 


(7) 


with the variation of t. Here, while calculating ai, 
we have subtracted £o which are ~ 2 and ~ 6.65, re¬ 
spectively, for Ti = oo and Tc. This subtraction is 
meaningful, considering the fact that the pure scaling 
with respect to time is contained m. £ — £q. Calculation 
of ai, without such subtraction, provided early time 
exponents much smaller than the theoretical expecta¬ 
tion for the conserved dynamics 0 . This has previ¬ 
ously been understood to be due to the curvature de¬ 
pendent correction. Such confusion has recently been 
corrected [l5|. Note that, there may be a delay time 
for a system to become unstable following a quench. 
Thus, for an appropriate understanding of a time de¬ 
pendent exponent, the value of £o need not be the 
length at t = 0. Via finite-size scaling analysis, this 
was demonstrated in a recent work |15j |. Here, how¬ 
ever, we do not undertake such a task. 

It may be argued that the value of £o should be of 
the order of system size for Ti = Tc, since ^ ^ L ai Tc. 
Note here that at criticality fluctuations exist at all 
length scales. At Ti = Tc, small value of £o, compared 
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to L, is due to the fact that our calculations did not 
ignore such fluctuations. 

First important observation from the inset of Fig. [5] 
is that the value of approaches 1 /2 from the upper 
side. This is in contradiction with the corresponding 
behavior for the conserved order parameter dynamics 
with Tf very close to zero 38|. In the latter case, 


the early time dynamics provides a growth exponent 
much smaller than the expected asymptotic value 1/3. 
Second, after t = 5, both the data sets practically 
follow each other, implying no difference in the growth 
of £ almost from the beginning! 



C being a constant. Such full forms are useful for a 
finite-size scaling analysis to accurately quantify the 
exponent a. It appears that even for a power-law 
behavior of f {!/£), the asymptotic behavior in the 
growth law is reached exponentially fast. Of course, 
from least square fitting exercise of the £ vs t data 
also one can aim to obtain the early time corrections. 
However, this method is more arbitrary. Often deriva¬ 
tives help guessing the functional forms better. We 
leave the exercise of a finite-size scaling analysis to 
accurately estimate a and the finite-size effects for a 
future work. 

Before moving on to presenting results in d = 3, 
we discuss the issue of persistence again. The essen¬ 
tial feature in the initial configurations prepared at 
different temperatures is the difference in the equilib¬ 
rium correlation length The basic question, prior 
to the study, one asks, how the size of ^ affects the 
decay of persistence probability? For each value of 
do we have a unique exponent 6 describing the full 
decay? The answer, as we have observed, is certainly 
not in affirmative. Essentially, the decay exponent for 
Tj = oo is recovered for all ^ (< oo) in the long time 
limit. The crossover to this asymptotic behavior gets 
delayed with the increase of It is then natural to 
ask if this crossover occurs when £ crosses This can 
be answered by appropriately estimating the crossover 
length £c. 


FIG. 5. Average domain sizes, £{t), are plotted vs t, for 
quenches to T/ = 0 from Ti = oo and Tc, in d = 2. The 
solid line represents behavior. Corresponding instan¬ 
taneous exponents, vs 1/t, are shown in the inset. All 
these results are for L = 512. 

From the length (or time) dependence of a^, one 
can write 


ai=a + f{i/£), ( 8 ) 


to obtain 


/ —f-—TTT = (9) 

J al[l + y{!/£)] 

If f{l/£) can be quantified accurately from the simu¬ 
lation data, a full time dependence of £ is obtainable. 
E.g., if f{l/£) is a power law, Apf£^, Ap being a con¬ 
stant, by taking al^ > Ap, one finds 


£l/a 1 


( 10 ) 


Assuming that a correction disappears fast, such that 
£ we obtain 


C 

al3t°^P 


), 


( 11 ) 



FIG. 6. Scaling plot of persistence probability versus t/tc 
where the crossover time (to the asymptotic behavior) tc 
has been used as an adjustable parameter to obtain opti¬ 
mum data collapse. Inset: Plots of — 1 versus e. The 
circles correspond to estimates of tc from G, the squares 
are directly obtained from the scaling plots of P vs ll£c- 
The solid line has d — 2 Ising critical divergence of corre¬ 
lation length. All results were obtained using L = 512 in 
d = 2. 


£ - exp( 
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In the main frame of Fig. [5] we show plots of per¬ 
sistence from different values of T^, for quenches to 
Tf = 0. Here the time axis is scaled by appropri¬ 
ate factors (expected to be proportional to cross over 
time tc) to obtain collapse of data in the asymptotic 
regime. Quality of collapse, on top of the Ti = oo data 
set, again confirms that 9 ~ 0.22 in the t ^ oo limit 
for all Ti (> Tc). From the square roots of these Ti 
dependent scaling factors, one can obtain £c which is 
expected to scale as £c ~ ^ ~ Note that for the 
Ising model v = 1 m. d = 2 and ~ 0.63 in d = 3. con¬ 
sidering that the Ti = oo data have been used as the 
reference case, it will be appropriate to fit the data 
set for ic to the form 

4 = l + dl,e-", (12) 

since ^ 1 for Ti = oo. Unless we are very close to 
Tc such additional term cannot be neglected. In the 
inset of Fig. [5] we ahve plotted £c — 1 as a function 
of e, on a log-log scale. The data set (circles) appear 
consistent with u = 1. When £c is extracted from 
tc, a better exercise requires incorporation of io and 
growth amplitude for each Ti. To avoid this problem, 
we have also obtained ic directly from the scaling plots 
of persistence data vs ijic (see the squares). Both 
data sets appear nicely consistent with each other. In 
fact, fitting these data to the form in Eq. © we 
obtain v ~ 0.95. 


B. d = 3 



FIG. 7. Plot of P vs t, for quenches from Ti = oo and 
Ti = Tc, to Tf = 0. In each of the cases results from 
two different system sizes are included. The solid lines 
have power-law decays with exponents 0.1 and 0.18, as 
indicated on the figure. All results corresponds to d = 3. 



FIG. 8. Instantaneous exponents 6i are plotted vs 1/i, 
for quenches from T = Tc to T/ = 0 in d = 3. Results 
from different values of L are included. The horizontal 
solid lines are related to the estimation of L-dependent 9. 
Inset: System size dependent 9 are plotted as a function 
of 1/L. The continuous line there is a linear fitting (see 
text for details). 


Next we present results in d = 3. All the impor¬ 
tant facts being discussed in the previous subsection, 
here we straightway present the results. Noting that 
nothing remarkable happened for domain growth in 
the lower dimension, we stick only to the persistence 
probability. 

In Fig. [3 we show the P vs t plots for quenches 
from Ti = oo and Ti = Tc, keeping Ty = 0 in both 
the cases. For each value of Ti, results from multiple 
system sizes are presented. The data for Ti = oo are 
consistent with 9 = 0.18, previously reported in 2l|. 
Thus, here we aim to accurately quantify the value of 
9 for T = Tc. 


Even though, for Ti = Tc, data from both the sys¬ 
tem sizes in Eig. [3 look consistent with each other, 
finite-size effects are detectable from a closer look. In 
the main frame of Fig. [3 we plot 9i versus l/£ for a 
few different values of L. Like in d = 2, from the flat 
regions we identify system size dependent 9, a plot of 
which is shown in the inset of this figure. Again, the 
9 vs 1/L data exhibits a reasonable linear trend and 
an extrapolation to L = oo provides 9 ~ 0.106. 

Similar to d = 2, for Tc < T; < oo, two step decays 
exist in d = 3 as well. In the main frame of Fig. 
[HKa) we have demonstrated the estimation of 9 for the 
first step, for two representative values of Ti. In the 
inset of Fig. [OKa) we have plotted these exponents 
as a function of e. A fit of this data set to the form 
in Eq. ([6|) provides 9{Ti = Tc) = 0.096, A = 0.07 
and X = 0.34. This value of 9 is in good agreement 
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this, as well as to get rid of the influence of Ti depen¬ 
dent io and growth amplitude, we have obtained ic 
from these plots only. 

Finally in the inset of Fig. IHJb) we plot £c — 1 as a 
function of e, on a log-log scale. The divergence of the 
length scale is consistent with a power law exponent 
0.63 which is the critical exponent for ^ in d = 3. In 
fact, a fitting to the form in Eq. (USD, excluding the 
point closest to Tc, provides v ~ 0.67. The reason for 
the exclusion of the last point is finite-size effect which 
is clearly visible. Data below e = 0.1 are, however, 
included in the nonlinear fitting in the inset of Fig. 
ini^a). This is the reason for quoting larger error bar in 
the final value of 6 there. 

IV CONCLUSIONS 



FIG. 9. (a) Estimation of 6 for the first step of decay is 
demonstrated in d = 3. Inset: Exponent 6 for the first step 
of the decay is plotted as a function of e, in d = 3. The 
continuous line is a non-linear fitting. Further details are 
provided in the text, (b) Scaling plot of P, versus tjic, in 
d = 3. The solid line has a power-law decay with exponent 
0.54. Inset: Plot of — 1, in d = 3, versus e. The solid 
line there has d = 3 Ising critical divergence. All results 
were obtained using L = 128. 


with the one obtained from Fig. |51 Thus, in d = 3, 
for Ti = Tc^ we quote 0 = 0.10 ± 0.02. The effect of 
growing correlation length in the initial configurations 
certainly appears weaker in this space dimension. 

In Fig. [9l[b) we show scaling plots of P, vs i/£c, us¬ 
ing data from different values of Ti. Collapse of data 
is good. The late time behavior is power-law and is 
consistent with a decay exponent 0.54. Considering 
that 0 cs 0.18 in d = 3, this implies a = l/3ind = 3 
(^ . As stated in Ref. [39], deviation of a, in this 
dimension, from 1/2, is not yet understood. To avoid 


In conclusion, we have studied phase ordering dy¬ 
namics in Ising ferromagnets for various combinations 
of initial (Ti) and final (Ty) temperatures in d = 2 
and 3. In this work, the primary focus has been on 
the persistence probability, P, and its connection with 
the growth of average domain size, £, as well as the 
equilibrium initial correlation length 

Our general observation has been that, irrespective 
of the value of Ti , the decay of P becomes faster with 
the increase of T/, after a certain critical number for 
the latter. This is understood to be due to spins af¬ 
fected by thermal fluctuations. When this effect is 
taken care of [l^, the long time decay appears to be 
power law with exponent 2^ 21 j consistent with the 
one for quench to P/ = 0. 

As Ti approaches Tc, two-step power-law decay be¬ 
comes prominent, the second part having exponent 
6 ~ 0.22 in d = 2 and ~ 0.18 in d = 3, same as 
Ti = oo and T/ = 0 case. For Ti = T^ thought to 
provide a new universality class, the first part of the 
two-step process lives for ever. The corresponding val¬ 
ues of the exponent have been identified to be ~ 0.035 
in d = 2 and ~ 0.10 in d = 3. Thus, even though, the 
decay of persistence probability is disconnected with 
the growth of domain length, its behavior is strongly 
connected with the initial correlation length. It has 
been shown that the crossover length to the second 
step of decays diverges as the equilibrium correlation 
length in both the dimensions. This leads to the ques¬ 
tion of difference in the fractal dimensions in the pre 
and post crossover regimes which needs to be looked 
at carefully. 

We have not observed any initial configuration de¬ 
pendence of the growth of the average domain size. 
This is consistent with a previous study 3^ but more 
explicitly demonstrated here. Essentially, even the 
transients are not affected due to change in initial 
temperature. However, stronger finite-size effects are 
detected for lower values of 7/. For domain growth. 




















a striking observation is that the early time exponent 
is much higher than the asymptotic value, despite Tf 
being zero. This is at variance with the conserved or¬ 
der parameter dynamics. These are all interesting new 
results, requiring appropriate theoretical attention. 

In future we will focus on persistence for the con¬ 
served order parameter dynamics. For the conserved 
dynamics, initial temperature dependence of aging 
and domain growth are also important open problems. 


ACKNOWLEDGEMENT 

The authors are thankful to the Department of Sci¬ 
ence and Technology, India, for financial support via 
grant number SR/S2/RJN-13/2009. 

* das@jncasr.ac.in 


[1] ONUKI A., Phase Transition Dynamics, (Cambridge 
University Press, Cambridge, 2002). 

[2] PURI S. and WADHAWAN V. , Kinetics of Phase 
Transitions, (Boca Raton, CRC Press, 2009). 

[3] BRAY A. J., Adv. Phys. 51, 481 (2002). 

[4] BINDER K., in Phase transformation of materials, 
edited by CAHN R. W., HAASEN P., KRAMER 
E. J., Mater. Sci. Technol. Vol. 5 (VCH, Wein- 
heim,1991), p. 405. 

[5] CROSS M. C. and HOHENBERG P. C., Rev. Mod. 
Phys. 65, 851 (1993). 

[6] OHTA T., JASNOW D. and KAWASAKI K., Phys. 
Rev. Lett. 49, 1223 (1982). 

[7] BRAY A. J. and PURI S., Phys. Rev. Lett. 67, 2670 
(1991). 

[8] TOYOKI H., Phys. Rev. B 45, 1965 (1992). 

[9] DAS S. K., PURI S. and CROSS M. C., Phys. Rev. 
E 64, 46206 (2001). 

[10] LIFSHITZ I. M. and SLYOZOV V. V., J. Phys. 
Chem. Solids 19, 35 (1961). 

[11] BINDER K. and STAUFFER D., Phys. Rev. Lett. 33, 
1006 (1974). 

[12] ALLEN S. M. and CAHN J. W., Acta Metall. 27, 
1085 (1979). 

[13] SIGGIA E. D., Phys. Rev. A 20, 595 (1979). 

[14] HUSE D. A., Phys. Rev. B 34, 7845 (1986). 

[15] MAJUMDER S. and DAS S. K., Phys. Rev. E 81, 
050102 (2010). 


[16] MAJUMDAR S. N., SIRE C., BRAY A. J. and COR¬ 
NELL S. J., Phys. Rev. Lett. 77, 2867 (1996). 

[17] MAJUMDAR S. N., BRAY A. J., CORNELL S. J. 
and SIRE C., Phys. Rev. Lett. 77, 3704 (1996). 

[18] DERRIDA B., HAKIM V. and ZEITAK R., Phys. 
Rev. Lett. 77, 2871 (1996). 

[19] DERRIDA B., Phys. Rev. E 55, 3705 (1997). 

[20] MANOJ G. and RAY P., Phys. Rev. E 62, 7755 

( 2000 ). 

[21] STAUFFER D., Int. J. Mod. Phys. C 8, 361 (1997). 

[22] CHAKRABORTY D. and BHATTACHARJEE J. K., 
Phys. Rev. E 76, 031117 (2007). 

[23] SAHARAY M. and SEN P., Physica A 318, 243 
(2003). 

[24] FISHER D. S. and HUSE D. A., Phys. Rev. B 38, 
373 (1988). 

[25] LIU F. and MAZENKO G. F., Phys. Rev. B 44, 9185 
(1991). 

[26] GORBERI F., LIPPIELLO E., ZANNETTI M., 
Phys. Rev. E 74, 041106 (2006). 

[27] AHMAD S., GORBERI F., DAS S. K., LIPPIELLO 
E., PURI S. and ZANNETTI M., Phys. Rev. E 86, 
061129 (2012). 

[28] MAJUMDER S. and DAS S. K., Phys. Rev. Lett. Ill, 
055503 (2013). 

[29] MIDYA J., MAJUMDER S. and DAS S. K., J. Phys. 
: Condens. Matter. 26, 452202 (2014). 

[30] DASGUPTA G. and PANDIT R., Phys. Rev. B 33, 
4752 (1986). 

[31] HUMAYUN K. and BRAY A. J., J. Phys. A : Math. 
Gen. 24, 1915 (1991). 

[32] BRAY A. J., HUMAYUN K. and NEWMAN T. J., 
Phys. Rev. B 43, 3699 (1991). 

[33] BLANGHARD T., CUGLIANDOLO L. F. and 
PIGCO M., J. Stat. Mech: Theory and experiment 
P12021 (2014). This very recent work was brought 
to our notice by an anonymous referee. The value 
6 = 0.035 obtained by us is extremely close to the 
conclusion from this study [3^ for triangular lattice. 
This implies, the lattice structure plays insignificant 
role. 

[34] SICILIA A., ARENZON J. A., BRAY A. J. and 
CUGLIANDOLO L. F., Phys. Rev. E 76, 061116 
(2007). 

[35] LANDAU D. P. and BINDER K., A Guide to Monte 
Garlo Simulations in Statistical Physics (Gambridge 
University Press, Gambridge, 2009). 

[36] GLAUBER R. J., J. Math. Phys. 4, 294 (1963). 

[37] FISHER M. E., in Gritical Phenomena, edited by 
GREEN M. S. (Academic Press, London, 1971). 

[38] MAJUMDER S. and DAS S. K., Phys. Ghem. Chem. 
Phys. 15, 13209 (2013). 

[39] GUEILLE S. and SIRE G., J. Phys. A: Math. Gen. 
30, L791 (1997). 



